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ABSTRACT 



Compared to their centimeter-wavelength counterparts, millimeter recombination lines (RLs) are intrinsically brighter and are free 
of pressure broadening. We report observations of RLs (H30o- at 231.9 GHz, H53or at 42.9 GHz) and the millimeter and centime- 
ter continuum toward the Becklin-Neugebauer (BN) object in Orion, obtained from the Atacama Large Millimeter/submillimeter 
Array (ALMA) Science Verification archive and the Very Large Array (VLA). The RL emission appears to be arising from the 
slowly-moving, dense (N e = 8.4 x 10 6 cnT 3 ) base of the ionized envelope around BN. This ionized gas has a relatively low elec- 
tron temperature (T e < 4900 K) and small (« 10 km s _1 ) bulk motions. Comparing our continuum measurements with previous 
(non)detections, it is possible that BN has large flux variations in the millimeter. However, dedicated observations with a uniform 
setup are needed to confirm this. From the H30cf line, the central line-of-sight LSR velocity of BN is 26.3 km s" 1 . 
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Understanding the nature of the radio sources in the Orion 
BN/KL region is important because these sources exhibit sig- 
nificant diverging proper motions that may trace the recent (~ 
500 yr ag o) disintegration of a compact cluster of young mas- 
sive sta rs (iRodrfguez et al.l2005l:lG6mez et al.l2008l:lGoddi et all 
I201 lbl) . This past interaction may also have caused the ex- 
plosive outflow seen in vibrationally-excited H2 and C O lines 
(IZapata et al.ll2009HBallv et al.ll201 ltlNissen et al.ll2012h . 

In particular, observing hydrogen recombination lines (RLs) 
can be used to obtain the line-of-sight velocities of the radio 
sources, which are needed to complement the proper-motion in- 
formation. Millimeter (mm) RLs are in principle better suited to 
giving the centtal velocity of the system because they are more 
optically thin, much less pressure broa dened, and intrinsically 
brigh t er than centimeter (cm) lines (e.g jGordon & Sorochenkol 
2002: IPeiters et al.ll2.012b . However, longer wavelength lines are 
a valuable complement because they are sensitive to the gas den- 
sity. By observing multiple RL transitions, hints to the structure 
of the ion ized nebula can be obtained even in unreso lved obser- 
vations ( Ke to et al . 2008: iGalvan-Madrid et al.ll2009l 

In this paper we report our findings in the recently released 
Science Verification (SV) data taken toward Orion BN/KL with 
the Atacama Large Millimeter/submillimeter Array (ALMA). 
We also report on previous RL observations obtained with the 
Very Large Array (VLA), as well as on the continuum obtained 
with both interferometers. 



2. Data 

Orion BN/KL was observed on January 19, 2012 with ALMA as 
part of the SV program. In this paper we make use of the fully 
calibrated visibilities released on April 12, 2012. The data con- 



sisted of 20 different spectral windows covering the frequency 
range from 213.72 GHz to 246.62 GHz. Each spectral window 
had a uniform spectral resolution of 0.488 MHz (~ 0.6 km s _1 ) 
and was observed for an on-source time of ~ 15 minutes. The 
array used 16 12-m antennas in a compact configuration with 
baselines in the range of 15 k/l to 269 kA. The pointing center 
was a J2 ooo = 05 h 35 m 14?35, feooo = -05° 22' 35^0. 

The absolute amplitude scale was derived from observa- 
tions of Callisto and Mars. The phase calibrator was J0607- 
085. We imaged the phase calibrator and measured a flux den- 
sity of S = 1 .40 Jy. The SMA data baseQ reports flux densi- 
ties for J0607-085 within ten days of the ALMA observation: 
So9Jan20i2 = 1.25 ± 0.06 Jy, 5 2 7Jan20i2 = 1.37 ± 0.07 Jy. From 
this, we estimate that the amplitude scale of the ALMA observa- 
tions is accurate to ~ 10 %. The processing was done in CAS A 
versions 3.3.0 and 3.4.0. 

We also present archival Very Large ArrajQ (VLA) data. The 
H53o' (vo = 42.95 19 7 GHz) B-array observa tions have been pre- 
viously presented by Rodrig uez et al.l (2009) (project AR635). 

Finally, several epochs of continuum observations at 8.4 
GHz in the A configuration were processed using standard meth- 
ods in the AIPS software. In all epochs 1331+305 was used as 
amplitude calibrator with an adopted flux density of 5.21 Jy. The 
epochs, project codes, phase calibrators, and their bootstrapped 
flux densities are given in Table Q] 

3. Results 

Continuum subtraction was done in the u-v domain by carefully 
selecting the line-free channels in each spectral window and fit- 



1 http://smal .sma.hawaii.edu/callist/callist.html/ 

2 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc. under cooperative agreement with the 
National Science Foundation. 
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Table 1. VLA calibrators at 8.4 GHz 



Epoch 


Project 


Phase 


Bootstrapped flux 






calibrator 


density [Jy] 


06 Sep 1991 


AM335 


0530+135 


1.77 ±0.04 


29 Apr 1994 


AM442 


0501-019 


2.38 ±0.01 


22 Jul 1995 


AM494 


0541-056 


1.97 ±0.14 


21 Nov 1996 


AM543 


0501-019 


1.51 ±0.01 


11 Jan 1997 


AM543 


0501-019 


1.45 ±0.01 


13 Nov 2000 


AM668 


0541-056 


1.11 ±0.01 


12 May 2006 


AR593 


0541-056 


1.39 ±0.01 



ting a linear spectral baseline. Continuum images from different 
spectral windows are consistent with each other, indicating a low 
contamination from unaccounted line emission. 

All the spectral windows were imaged with a threefold pur- 
pose: to detect the H30ff RL in the BN/KL radio sources (I, BN, 
and n), to detect molecular emission from BN, and to make a 
continuum image. Whereas we succeeded in detecting the RL 
in BN, it was not detected elsewhere. Also, no molecular emis- 
sion was detected from BN. We present maps that have not been 
corrected for primary beam attenuation. The corrected fluxes are 
quoted in the text and tables. 

3.1. Continuum. 

The low declination of Orion BN/KL, as well as the bright, ex- 
tended emission dominated by the "hot core" (HC), limit the dy- 
namic range of synthesis imaging and complicate the detection 
of weak sources in this field at (sub)mm wavelengths. Figure [TJ 
shows a deep continuum image from the ALMA data obtained 
by co-adding the line-free continuum visibilities from all spec- 
tral windows. The average frequency (wavelength) in this map 
is 230.17 GHz (1.30 mm). The dynamic range in this image is 
determined by limitations in the u-v coverage, and it is in the 
range of ~ 10 to 15 mJy. Our continuum map is ~ x2 deeper 
than the released continuum imag fl and similar in quality to the 
CARMA 1.3 -mm map reported b ylFriedel & Snvderld2008l). 

Th e centimeter radio source I dReid et al.ll2007t iGoddi et"aD 
1201 lbl) is only ~ 0"7 from the 1.3-mm peak of the HC, which 
has an extension of ~ 8"6 x 3"2 (see Fig. [TJ. Therefore, the 
continuum emission of source I cannot b e resolved from the 
HC in arcsecond-r esolution images (e.g.. iBeuther et ail [2005: 
IGoddi et"aTI 1201 lab . The BN object, as traced by its cm con- 
tinuum, is ~ 6"7 from the hot-core continuum peak. In the map 
shown in Fig. [TJ we detect BN with a peak intensity of 55 mJy 
beam -1 (112 mJy beam -1 when primary beam attenuation is 
taken into account). However, it is still at the intensity level of 
some of the extended emission. 

To disentangle the compact continuum emission from the ex- 
tended emission in the field, we made cleaned maps after re- 
moving the short baselines. This degradation of the u-v cover- 
age produces artifacts (closely-spaced sidelobes) close to source 
I and the HC, and we refrained from using this part of the image. 
However, at the position of BN, where emission is fainter, the 
restricted u-v coverage effectively removes the extended emis- 
sion around BN. We found that a minimum baseline length of 
as 100 k/l gives optimal results, filtering out the emission on 
scales > 1"4. Figure [2] shows this continuum map around BN. 
The noise is x5 lower than in Fig. 1 (~ 3 mJy). The 1.3-mm 
flux density of BN in this image is 62 mJy. After primary beam 
correction, the flux density of BN is S i.3 mm = 126 mJy. 

3 http://almascience.eso.org/alma-data/science-verification 
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Fig. 1. Continuum ALMA image centered at 230.17 GHz (1.3 
mm) of the Orion BN/KL region. The centimeter sources I, BN, 
and n are marked with circles. The emission is dominated by 
the bright, extended emission from the Orion "Hot Core", with 
a peak intensity of 834 mJy beam -1 . The rms noise is <x ~ 10 to 
15 mJy. The half-power beam width (HPBW) is i:'49 X V.'Ol, 
with a position angle PA= -3.2°. 

We note that the flux density for BN is the same, within 20 %, 
regardless of the u-v coverage restriction applied. A faint com- 
pact source with primary beam corrected flux S i.3 mm = 26 mJy 
is also detected 4"3 south of BN. This newly detected source is 
possibly ass ociated with the exten ded mid-infrared sources IRc6 
and IRc6 N (Shupi ng etal.ll200l . 
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Fig. 2. Continuum ALMA image centered at 231.17 GHz (1.3- 
mm) of the compact emission of the BN object. Only baselines 
with length B > 100 kA were used to filter out the extended 
emission in the field. The BN object is detected as an unresolved 
source. The rms noise around BN in this map is <x ~ 3 mJy. 
HPBW= 1!'32 x 0:'62, PA= 175.5°. 

It is possible that BN is variable in the (sub)mm contin- 
uum. While our detection of BN is consistent with the 218 GHz 
flux of S = 150 + 50 mJy (from OVR O observatio ns in 1993 
and 1994) reported bv iPlambeck et alJ (1 19951) and iBlake et all 
(fl9%|) . lFriedel & Snvderi d2008l) do not detect BN at 231 GHz 
down to a 3cr level of ~ 34 mJy (CARMA observations in 2007). 
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Table 2. BN fluxes. 



0.20 



Epoch 



[mJy] S 1.3mm [mJy] Notes 



06 Sep 1991 
Nov 1993 to Feb 1994 
29 Apr 1994 
22 Jul 1995 
21 Nov 1996 

11 Jan 1997 
13 Nov 2000 

12 May 2006 
Mar 2007 

19 Jan 2012 



4.61+0.20 

3.73 ± 0.04 
3.36 + 0.3 
3.43 + 0.1 
3.90 ± 0.07 
3.83 + 0.14 
4.78 + 0.11 



150 + 50 



< 34 
126+14 



Notes, a: VLA- A archival data, HPBWk 0"3 (this paper), b: 
OVR O observations at 21 8 GHz, HPBW= V.'S x V.'O (Plambec keTaTl 
19951 : iBlake et alj 1 19961) . c: VLA-A, HPBW« 0'.'3 dZapata et alJ 
20041). d: CARMA obs ervations at 231 GHz, HPBW= 2'.' 5 x 01' 8 
dFriedel & S nyder 200|). e: ALMA SV data, HPBW= L'3 x 0'.'6 (this 



paper). The error is determined by variations in the obtained flux with 
different apertures in the S rain = 100 kA map. The flux density has been 
corrected for primary beam attenuation. 



However, the significance of these variations can be affected 
by the inhomogeneity of the mm observations. We have repro- 
cessed three epochs of VLA archival data at 3.6 cm (8.4 GHz) to 
look for possible variations in the continuum emission from BN. 
Table |2] lists the me asured 3.6-cm fluxe s, together with the four 
epochs reported by Zap ata et al] (|2004), also at 3.6 cm, and the 
available 1-mm band measurements. Based on this data set, the 
3.6-cm flux is quite stable (within 30 % of the maximum value). 

3.2. Recombination lines 

We report the first detection of mm RL emission toward the 
BN object. We also searched for H30a emission in the rest 
of the Orion BN/KL field, in particular from sources I and n, 
without success. Just as with the continuum, the proximity of 
source I to the HC contaminates any emission from the rel- 
atively faint RL with brighter, more extended molecular lines 
with rest frequencies within a few MHz of the H30q- line. Figure 
3 shows the H30a and the H53a (data originally presented in 
iRodrfguez et al.l 120091) spectra toward BN. The line parame- 
ters are summarized in Table [3] Figure 4 shows the velocity- 
integrated intensity maps (moment 0) for both lines. Similar 
to the mm continuum, the H30a map uses only baselines with 
length > 100 



Table 3. Recombination line parameters. 



Transition 


y rest 


S peak 


VLSR,peak 


FWHM 




[GHz] 


[mJy] 


[kms- 1 ] 


[km s" 1 ] 


H30a 


231.9009 


181 ± 14 


26.3 ± 0.5 


15.4+ 1.3 


H53a 


42.95197 


9.9+1.1 


19.8 + 2.1 


39.0 ± 5.0 



Notes. From Gaussian fits to the spectra shown in Fig. 3 using an aper- 
ture adjusted to the emitting area shown in Fig. 4. 

The roOff line is narrower than the HSSff line, and it is red- 
shifted by a few km s _1 . This behavior has been observed in hy- 



4 The H30o- line is also detected in maps that use the complete base- 
line range, but bright molecular emission from the HC centered at ~ 10 
MHz blueshifted with respect to the RL creates sidelobes that are for 
the most part removed by taking only the longest baselines. 




20 40 

Vlsr [km/s] 

Fig. 3. Hydrogen recombination line spectra (gray) and 
Gaussian fits (red) toward BN. The top panel shows the H30a 
line observed with ALMA. The flux is corrected by primary 
beam attenuation. The bottom panel shows the H53a line ob- 
served with the VLA. The channel widths are 2.0 km s _1 and 
10.9 km s _1 , respectively. 

percompact HII regions around embedded y oung massive stars 
ilKeto et alJl2008k iGalvan-Madrid et al.ll2009h . Although BN is 
not deeply embedded in molecular gas, the same interpretation 
seems to apply: that the difference in linewidth is caused by the 
action of "pressure" broadening from particle collisions in the 
low-frequency line, whereas the velocity differences are caused 
by optical depth effects. In an outflowing gas where the denser 
part is closer to the central star, the optically thinner, higher fre- 
quency line will trace the centermost gas from both the front and 
the back sides of the nebula. In contrast, the optically thicker, 
lower frequency line, will trace gas farther away from the cen- 
tral s tar, mainly from the front (blueshifted) side of the nebula 
(see lKetoetai] f2008). Therefore, the H30o' line should give a 
better estimation of the line-of-sight velocity of the star inside 
the BN nebula, v L sr(BN) = 26.3 + 0.5 km s" 1 . 

The total RL width (Av) has contributions from thermal 
broadening (Avth), dynamical broadening from macroscopic gas 
motions (Avay), and pressure broadening (Av pr ). We elaborate 
this aspect in Appendix A. The line FWHM of the combin ed 
profile is given by (equation 2 of IGalvan-Madrid et al . 2009): 



Av 



0.534Av pr + (Av 2 dy + Avl + 0.217Avj; r ) 1/2 . 



(1) 



The H30q- line is only broadened thermally and by dynam- 
ical motions (Av pr < 0.6 km s _1 from equation A. 4). Also, 
in the BN object this line appears to be too narrow (15.4 km 
s -1 ) even for pure thermal broadening at T e = 10 4 K (equation 
A. 2). Therefore, for the gas traced by the H30a line, T e is at 
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J2000 Right Ascension 

Fig. 4. Maps of the velocity-integrated (moment 0) emission of 
the recombination lines. The color scale shows the LOOa line 
observed with ALMA (HPBW= L'40 x 0!'63, PA= 177.9°). 
Contours show the H53a line observed with the VLA (HPBW= 
0:'22xO"18,PA= 19.4°). 

most 4900 K, and the bulk gas motions are constrained to be 
m uch smaller than the thermal width. Finally, using equation 3 
of lKetoetai] (12008). the average electron density traced by these 
RLsis/V c = 8.4 X 10 6 cm 3 . 

Non-LTE effects and dust opacity may play a role for 
(sub)mm RLs. In Appendix B we show that the RLs in BN are 
reasonably close to LTE and that it is possible that dust in the 
line-of-sight makes a considerable contribution to the 1.3-mm 
continuum. 

4. Discussion 

One possible explanation for the differences in the mm fluxes 
that have been reported for BN is that they are, at least partially, 
caused by systematic errors in the measurements. All the com- 
piled cm fluxes are measured with the VLA A-array, which has 
an absolute flux-scale uncertainty of a few percent, while the 
flux scale of mm interferometers is known to be only accurate 
to within -10 9^ Differences in u-v coverage may also mimic 
flux variations, but since the continuum emission from BN is 
less than one arcsecond in size, this effect should not be im- 
portant. Indeed, the flux that we recover from BN is the same 
within ~ 20 % for different u-v restrictions, (see Section 3.1). 
Alternatively, if BN is variable in the mm, it could be caused by 
the sudden recombination of dense gas (the timescale of whic h 
is ~ 1 month for densities N e ~ 10 cm , Oster brockiri989h . 
We emphasize that only dedicated multiepoch observations with 
a uniform setup can clarify this issue. 

The BN object is known to b e a young star with a mass M* » 
10 M (Rodrig uez et~aTll2005l) . In the absence of other forces 
(e.g., wind acceleration mechanisms), this mass is enough to 
gravitationally confine ionized hydrogen within a radius of « 55 
AU (see lKetoll2007l) . If the H30a- is emitted within a volume as 
small as the radio continuum (radius ~ 25 AU, Rodriguez et al.1 
2009), this could explain the lack of bulk motions in this RL. 
The H53ff line appears to be tracing gas that is outflowing to- 
ward us with a velocity ~ 6.5 km s _1 blueshifted with respect to 

3 http://www.vla.nrao.edu/astro/calib/manual/, 
http://haneul.astro.illinois.edu/~wkwon/CARMA/fluxcal/ 



the systemic. Infrared and optical RLs have also been detected 
towar d BN and extend to velocities > 100 km s _1 (IScoville et al.l 
1983). However, these RLs are significantly affected by obscura- 
tion and non-LTE effects, so they may trace gas from a different 
volume than the radio and (sub)mm RLs. 

5. Conclusions 

We report the first detection of 1.3 mm RL (LOOa) emission to- 
ward Orion from the BN object. The mm RL appears to be aris- 
ing from the dense (N e = 8.4 x 10 6 cm -3 ), static (bulk motions 
<< 10 km s _1 ) base of the ionized nebula around the central 
massive star in BN, whose line-of-sight LSR velocity is now es- 
timated to be vlsr(BN) = 26.3 + 0.5 km s -1 . Compared to the 
H30or line, the optically thicker H53o' seems to trace gas out- 
flowing toward us at ~ 6.5 km s _1 . 

Dedicated ALMA observations at subarcsecond resolution 
would be able to detect RL emission from the rest of the Orion 
BN/KL sources, to obtain their 3D velocities by combining them 
with the available proper motion measurements, and to quantify 
their possible flux variations. 
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Appendix A: Linewidth of recombination lines 



Appendix B: cm and (sub)mm recombination lines 



RLs are broadened by several mechanisms: a natural broaden- 
ing from the quantum uncertainty AE of the energy level E, 
a thermal/microturbulence Gaussian broadening from the mo- 
tions of the emitting particles and of parcels of gas that are much 
smaller than the beam, the Stark (or "pressure") broadening from 
the perturbation of the atomic energy levels by the electric field 
of neighboring charged particles, and a dynamical broadening 
from bulk flows (e.g., infall, rotation, outflow) in the gas. A de- 
tailed discu ssion of these processes and t he physics of RLs can 
be found in lGordon & Sorochenkol(l2002l) . 

For n > 20 (where Hna is a n + 1 — > n transition), the natural 
linewidth in velocity units is 



Av N ,1.2xl0V n( " + 1) 



(n+1) 



2 ' 



(A.l) 



where c is the speed of light. For H30or and HSSa, Avn is 1.3 x 
10~ 3 km s _1 and 0.5 x 10~ 3 km s _1 , respectively. The natural 
broadening is negligible for radio and (sub)mm RLs. 

The thermal distributions of velocities of both the individual 
particles and small pockets of gas ("microturbulence") produce 
a Gaussian contribution to the broadening. Neglecting microtur- 
bulence, the thermal FWHM is 



Av f 



81n2 



m H 



1/2 



(A.2) 



where is the Boltzmann constant, is the hydrogen-atom 
mass, and it is assumed that all the gas is thermalized to the elec- 
tron temperature T e . For a pure-hydrogen gas with T e = 9000 K, 
Av th = 20.3 km s" 1 . 

The pressure broadening has a Lorentzian shape, and it in- 
creases with density and quantum number. Collisions with ions 
and electrons make differen t contributions. The c ontribution 
from ions has the form dGordon & Sorochenk o 2002) 

Ar n „i - | /V, — 1(0.06 + 2.5 x lO^rjI-r^-) fl + ^^\(A.3) 



, p _k.£(0.06 + 2.5xlO-T e )^ 



n + 1 



where y, = 6 - 2.7 x lO^ - 0.13(« + 1)/100. For T e = 9000 
K and N, = 10 7 crrT 3 , the H30a line is virtually free from ion 
broadening (0.04 km s _1 , and decreases close to linearly with 
decreasing T e ), whereas the HSSff line is broadened by 5.1 km 
s" 1 . 

The collisions with electrons dominate the collisions with 
ions under most cond itions. They produce a broadening of width 
(ISmirnov et al.ll 1984 : 



Av, 



pr,e 



8.2JV, 



n + 1 
v A 100 



4.5 



2.25 An\ 
n + 1 }' 



(A.4) 



For N e = 10 7 cm 3 , the electron broadenings for the H53or and 
H30q' are « 37.3 km s and 0.6 km s , respectively. 

Finally, the last source of broadening is bulk motions (Avjy) 
of the ionized gas, which may be in the form of outflows/winds, 
inf all/accretion, and/or rotation. For unresolved observations, 
the method outlined here permits estimating the magnitude of 
these motions from the nonthermal linewidth of the line that is 
free of pressure broadening. Of course, ideally one wishes to 
spatially resolve these motions. ALMA in its Cycle 1 is now 
able to resolve the ionized-gas motions in sources as faint as 
~ 10 mJy at the line peak at subarcsecond resolution. 

Including all the contributions to the broadening, the total 
linewidth will have a Voigt profile with FWHM given by equa- 
tion (1) in the text. 



Here we discuss the advantages and limitations that cm and 
(sub)mm RLs have compared to each other. The cm lines are 
intrinsically fainter and optically thicker, and can be partially 
absorbed by the r elatively high contin uum opacity at their wave- 
lengths (see e.g., [Wilson et aD l2009h . A possible drawback of 
(sub)mm RLs is that the free-free continuum may be con- 
taminated by dust emission. Another problem is that they are 
out of LTE more easily than cm li nes, so their interpretation 
may require care ful modeling (e.g.. Uimenez-S erra et alj|201 U 
iPeters et alj[2012l) . As shown below, our simple LTE interpreta- 
tion seems to be reasonable, but it is possible that a significant 
fraction of the 1.3-mm continuum comes from dust in the line of 
sight. 

Let b„ be the ratio of the population of level n to its LTE pop- 
ulation. Then the actual line abso rption coefficient /cl is relat ed 
to the LTE coefficient kl.lte by dGordon & Sorochenko 2002) 



with 

" b„ \ hv dn 



(B.l) 



(B.2) 



While < b n < 1 expresses the non-LTE level depopulation, /3„ 
can be negative when there are conditions for maser amplifica- 
tion. 

In the optically-thin case (Tc,i.3mm ~ 0.08 from the data), the 
non-L TE corrected line intensity /l is dGordon & Sorochenko! 
12002b 



/ L ~/L,LTE& n+ i(l-T C j3/2). 



(B.3) 



Us ing the tabulated values in the calculations of Wal mslevl 
d 19901) . b 30 * 0.98 and £30 ~ -1.5 for a density N e = 10 7 cm" 3 . 
Therefore, the intensity of the MOa line is only amplified by 
~ 3 %. The ^3^ is even closer to LTE, with b 5i * 0.999 and 

,053 ~ 0.5. 

Sub(mm) RLs are more opticaly thin and intrinsically 
brighter than cm RLs. For a given T e and N e , and in the 
Rayleigh- Jeans regime (hv « k^T e ), the line absorption co- 
efficient decreases linearly with frequency °c v _1 , so the LTE 
line emission coefficient j'l = KLB v (T e ) oc v. Similarly, the con- 
tinuum absorption coefficient kq oc v -21 , so the LTE continuum 
emissivity jc v~° 1 . Therefore, in the optically thin regime, the 
line-to-continuum ratio increases almost linearly with frequency 
S^Av/Sc v 11 . Under these assumptions, the expected ratio 
of the H30or to HSSo- line-to-continuum ratios is 6.4. However, 
the observed value is 1.6. If we assume that only ~ 32 mJy 
of the 1.3-mm continuum are due to free-free then the line-to- 
continuum ratio would increase to the expected value. This sug- 
gests that dust contributes a significant fraction to the (sub)mm 
continuum of BN. However, it is known (and we have confirmed 
this with the ALMA data set) that BN is not embedded in dense 
molecular gas. Therefore, it is unlikely that this line-of-sight dust 
arises in a dense core around BN. 



